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Abstract 

We present the results of optical panoramic and long-slit spectroscopy of the neb- 
ula MF16 associated with the Ultraluminous X-ray Source NGC6946 ULX-1. More 
than 20 new emission lines are identified in the spectra. Using characteristic line 
ratios we find the electron density n e ~ 600 cm -3 , electron temperature in the range 
from ~ 9 000 K to ~ 20 000 K (for different diagnostic lines) and the total emitting 
gas mass M ~ 900 M Q . We also estimate the interstellar extinction towards the 
nebula as Ay ~ l m . 54 somewhat higher than the Galactic absorption. Observed line 
luminosities and ratios appear to be inconsistent with excitation and ionization by 
shock waves so we propose the central object responsible for powering the nebula. 
We estimate the parameters of the ionizing source using photon number estimates 
and Cloudy modelling. Required EUV luminosity (~ 10 40 erg s _1 ) is high even if 
compared with the X-ray luminosity. We argue that independently of their physical 
nature ULXs are likely to be bright UV and EUV sources. It is shown that the UV 
flux expected in the GALEX spectral range (1000—3000 A) is quite reachable for UV 
photometry. Measuring the luminosities and spectral slopes in the UV range may 
help to distinguish between the two most popular ULX models. 

Key words: ISM: individual (MF16), jets and outflows, supernova remnants X- 
rays: individual (NGC 6946 ULX-1) ultraviolet: general 
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1. Introduction 



1.1. ULXs and their Optical Nebulae 

A point-like X-ray source in an external galaxy is considered an Ultraluminous X-ray 
source (ULX) if its luminosity exceeds 10 39 erg s _1 and it is not an active galactic nucleus. 
It also makes sense to exclude X-ray bright SNe and young (such as several years) Supernova 
Remnants (SNR) that can shine as bright as 10 41 erg s _1 in X-rays like the remnant of SN1988Z 
(Fabian & Terlevich 1996). There are more than 150 ULXs known at the present time (Swartz 
et al. 2004) but very little is clear yet about the physical nature of these objects. 

ULXs are widely accepted as a class of accreting compact objects violating the Eddington 
luminosity limit for a conventional stellar mass black hole (about 1.3 x 10 39 erg s _1 for 10 M Q ). 
High luminosity may be a consequence of a higher accretor mass, supercritical accretion, mild 
geometrical collimation, relativistic beaming (Kording et al. 2002) or some combination of these 
effects. Two models are the most popular: (i) Intermediate-Mass Black Holes (IMBHs) with 
masses in the range 10 2 — 10 4 M (Colbert & Miller 2005; Madau & Rees 2001) accreting at sub- 
Eddington rates; (ii) supercritical accretion discs like that of SS433 around stellar mass black 
holes observed at low inclinations i < 20° (Katz 1986; King et al. 2001; Fabrika & Mescheryakov 
2001). Some authors (Soria & Kuncic 2007) propose "hybrid" models involving 20 — 100 M 
black holes accreting in a moderately supercritical regime. Supercritical accretion discs are 
considered by different authors either in the wind-dominated regime first introduced by Shakura 
& Sunyaev (1973) or in the advection-dominated regime known as slim disc (Abramowicz et 
al. 1988; Okajima et al. 2007). 

Large number of ULXs are surrounded by large-scale bubble nebulae (Pakull & Mirioni 
2003), probably shock-powered. ULX nebulae (ULXNe), however, do not form a homogeneous 
class of objects. In Abolmasov et al. (2007) we review the observational properties of 8 ULXNe 
including the object under consideration. We conclude that only about 50% of ULXNe may 
be considered shock-powered shells. In some of the observed ULXNe high-excitation lines such 
as [O III] A4959, 5007 are enhanced. In some cases HeIIA4686 emission is detected as bright as 
0.2—0.3 in H/3 units (Lehmann et al. 2005). The line may have stellar origin (Kuntz et al. 2005) 
as well as nebular. Measuring line widths may help to distiguish between these two cases. The 
nebulae may clarify the nature of corresponding X-ray sources in two ways: (%) probing the 
photoionizing radiation of the central object via gas excitation and ionization conditions; (ii) 
detecting and measuring jet/wind activity by kinematical effects and shock emission. 

Dynamically disturbed gas may be considered an evidence against the IMBH hypothesis 
because a standard accretion disc is unlikely to produce strong jets or wind. Supercritical 
accretion, on the other hand, is supposed to provide a massive outflow in the form of a strong 
wind carrying practically all the accreted mass and having kinematical luminosity ~ 10 38 erg s" 1 
in the case of SS433 (Fabrika 2004). Power comparable with the Eddington luminosity (~ 
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10 39 erg s _1 ) may be provided in the form of relativistic jets. Their mechanical luminosity 
acting for ~ 10 5 years is sufficient to produce a wind-blown bubble with properties close to 
those of some of the existing ULX nebulae (Pakull & Mirioni 2003; Abolmasov et al. 2007). 

ULXNe exhibit very diverse observational properties such as size, morphology and line 
ratios (Abolmasov et al. 2007; Pakull & Mirioni 2003). Integral luminosities in Balmer lines are 
in some cases of the order 10 38 erg s" 1 or higher that requires an energy source with luminosity 
> ^q39 er g s -i i nc li ca ting that ULXs are really powerful objects. Therefore strong beaming 
effects are excluded as the reason of apparently high luminosities of ULXs. For some objects 
like HoIX X-l and IC342 X-l the nebula is clearly shock-powered. In other cases (HoII X-l, 
M101-P098) low intensity of low-excitation lines, bright [O III] and He II emissions and quiet 
kinematics 1 suggest photoionization to be the main energy source. 

The only proven example of a persistent supercritical accretor is the peculiar binary 
SS433 (Fabrika 2004). Its accretion rate is of the order 10~ 4 M Q yr _1 , indicating mass transfer 
on the thermal timescale of a massive star. Most of the accreting material is ejected in the 
form of optically-thick wind with velocities V ~ 1000 — 2000 km s _1 . The total kinetic power 
of the jets of SS433 is of the order of 10 39 erg s _1 . SS433 is surrounded by a large (roughly 
40 x 120 pc) radionebula W50 (Dubner et al. 1998). The nebula harbors several optical emission- 
line filaments situated in the propagation direction of the jets. Large difference in angular size 
and heavy absorption in the case of W50 make comparison with ULXNe difficult. 

1.2. NGC6946 ULX-1 

The X-ray source under consideration is known as NGC6946 ULX-1 (Swartz et al. 2004), 
NGC6946 X-8 (Lira et al. 2000), N58 (Holt et al. 2003) or NGC6946 X-ll (Roberts & Colbert 
2003). It was first detected by ROSAT (Schlegel 1994) and identified with an optical nebula 
by Blair & Fesen (1994) and a marginally resolved radiosource by van Dyk et al. (1994). 

The nebula is known as MF16 - Matonick and Fesen 16 (Matonick & Fesen 1997). It 
was considered a SNR luminous in X-rays, in fact the most luminous in X-rays among the 
optically bright SNRs (Dunne et al. 2000), till it was proved by Roberts & Colbert (2003) that 
the X-ray emission originates from a much more compact source in the center of the nebula. 
Chandra source coordinates for J2000 epoch are: a = 20 h 35"W.74, 5 = +60° 11'30"6. 

The host galaxy NGC6946 is a late-type spiral with active star formation (de Gioia- 
Eastwood & Grasdalen 1984) and the greatest number of detected supernovae (no less than 8, 
according to Li et al. (2005)). The distance to the galaxy is estimated by different authors as 5.1 
(de Vaucoulers 1978), 5.5 (Tully 1988), 5.7 (Eastman et al. 1996) and 5.5 Mpc (Karachentsev 
et al. 2000) therefore we assume D = 5.5 Mpc. The divergence in distance estimates leads to a 
15% uncertainty in all the luminosities inferred. Spatial scale for D = 5.5 Mpc is about 27 pc 

1 For HoII X-l Lehmann et al. (2005) report velocity dispersion ~ 13 km s _1 in the vicinity of the X-ray 
source. 
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per arcsecond. 

A point source in the center of MF16 was detected by Blair et al. (2001) in HST data. 
The source (star d as denoted by Blair et al. (2001)) has a V magnitude of 22 m . 64 and a colour 
index of B — V = II 146. The object is significantly redder than the nearby stars a-c offset to 
the North- West by several seconds. If one accounts only for Galactic interstellar extinction 
(see discussion in section 3.2) equal to A^y AV) = l m .14 according to Schlegel et al. (1998) in 
the direction of MF16, the point-like counterpart resembles an early A la supergiant with 
My ~ — 7 m and B — V ~ m . Additional intrinsic absorption Ay ~ 0^5 (reported by Blair et 
al. (2001)) places the object slightly above the la sequence in the upper left part of the HR 
diagram. 

Optical emission-line spectrum of the nebula is generally consistent with the suggestion of 
shock heating but also contains high-excitation lines such as He IIA4686 and [O III] A5007, 4959 
doublet (Blair & Fesen 1994; Blair et al. 2001) too intensive for optically-bright SNRs (Matonick 
et al. 1997). Partially radiative shocks were proposed to explain the emission-line spectrum 
of MF16 (Blair et al. 2001) but they require a very powerful energy source (see discussion in 
section 4). 

The object is much brighter in the optical than a usual SNR. Its Ha luminosity is 
about 2 x 10 38 erg s -1 , an order of magnitude higher than the upper limiting Ha luminosity 
for optically-bright SNRs in nearby galaxies (Braun & Walterbos 1993). Dunne et al. (2000) 
detected two-component structure in the emission lines of the object in high-resolution echelle 
spectra. Broader components suggest expansion velocities about 250 km s _1 , narrower compo- 
nents have velocity dispersion ~ 20 — 40 km s" 1 . 

MF16 is also known as a radio source as bright as 1.3mJy at 20cm (van Dyk et al. 
1994) marginally resolved by VLA observations as a ~ 1" size object. The radio source appears 
to be displaced relative to the X-ray object and the optical source d by ~ 0"5 (see figure 1). 
Luminosity of the radio counterpart is about 20 times higher than that of W50 at 20cm (Dubner 
et al. 1998). The size of the radiosource is poorly known but the optical nebula is about 5 times 
more compact than W50. 

MF16 is practically isolated from other HII regions. The nearest one with comparable 
brightness is situated about 200 pc away (Blair et al. 2001). The inner brighter part of the 
nebula has a shape of a shell 20 pc x 34 pc elongated in East to West direction with a brighter 
Western loop. In the deep HST images a faint asymmetric halo may be seen (Blair et al. 2001). 

In the following section we describe the observational data. In section 3 we present the 
results of spectral analysis. In section 4 we analyze the excitation and ionization sources of the 
nebula and model the observed spectrum using Cloudy. In section 5 we discuss our results and 
the perspectives of UV observations of ULXs. 
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2. Observations and Data Reduction 



All the data were obtained at the Russian SAO 6m telescope. The main information 
about the observational material is listed in table 1. 

The data were flat-fielded and wavelength-corrected. Wavelength calibration was made 
using a He-Ne-Ar lamp. Atmospheric extinction was corrected using the data obtained by 
Kartasheva & Chunakova (1978) fitted by a simple formula Am = (0.013(A/10 4 A) _4 + 0.13)secz. 

2.1. Integral-Field Spectroscopy 

The MultiPupil Fiber Spectrograph (MPFS) is described in Afanasiev et al. (2001). 
MPFS is a panoramic spectrograph with 16" x 16" field of view consisting of 256 1" x 1" spatial 
sampling elements. Spectra from individual elements arranged according to their coverage of 
the celestial sphere form a three-dimensional structure that is usually called a data cube. Extra 
17 fibers collect sky background from 4' away offset regions. 

Data reduction system was developed in IDL environment. The reduction process in- 
cludes the procedures standard for panoramic spectroscopy data reduction (see for example 
Sanchez (2006)): bias subtraction, cosmic hits removal, flat-fielding, individual fiber extraction, 
wavelength and flux calibrations. The dispersion curve has an accuracy of about 5 — 10 km s -1 . 
Wavelength shifts between different fibers were corrected using the brightest night sky lines 
[OI]AA5577,6300,6364. We used GD248 spectrophotometric standard (Oke 1990) for flux cali- 
bration. Sensitivity variations of individual fibers were corrected using twilight sky exposures. 

We also added atmospheric dispersion (or atmospheric differential refraction, ADR) 
correction with accuracy better than 0."1 in order to achieve correct relative astrometry for the 
datacube. We determined the value of ADR using the spectral standard star. Its shift along 
the vertical as a function of wavelength was approximated by a second-order polynomial that 
resulted in a roughly 0"1 scatter around the best-fit parabola and no detectable systematical 
deviations. The spectral standard was observed at different zenith distance (z s tar —40°, z Q bj — 
55°) but we corrected for that difference using the well-checked octanz dependence for ADR. 
According to Filippenko (1982) octane approximation holds within ~0.2% at the 2 km altitude 
in the optical range up to secz ~ 5 (z ~ 80deg). 

2.2. Long-Slit Spectroscopy 

The observations were held with SCORPIO focal reducer (Afanasiev & Moiseev 2005), 
in the long-slit mode with moderate spectral resolution and slit width 1" . The position angle 
[PA = 96°.5) of the slit was chosen along the major axis of the nebula (figure 1). Reduction 
was performed in IDL environment and includes all the standard reduction steps for long-slit 
spectroscopy. 

The signal-to-noise ratio for the long-slit data is about factor 2-3 higher than in the 
MPFS data, mainly because of the better spectrograph transparency. Integral field data, how- 
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Fig. 1. V-band image obtained during the observations with SCORPIO (on the left panel) and HST ACS 
archival image (F658N filter) of the source vicinity (right panel, the region is shown on the left panel 
by a dashed rectangle). MPFS field is shown by a square. Slit position is shown on the left panel by 
two horisontal lines. Point sources c and d are marked. White circle is centered at the position of the 
X-ray source, black contours correspond to the VLA radio isophotes at 4.86GHz (van Dyk et al. 1994). 

ever, have the advantage of giving flux values free from slit losses and radial velocities without 
systematical errors usual for long-slit spectroscopy. Besides this, MPFS data cover larger spec- 
tral range. 

3. Spectroscopy Results 

3.1. Emission Line Spectrum 

MPFS spectrum was extracted within 2" radius diaphragm, the long-slit spectrum 
was integrated within a ±l!'5 band. Emission line fluxes were calculated using gauss anal- 
ysis. Fitting with two gaussians was used for [SII]AA6717,6731 doublet, H7+[OIII]A4363, 
HeIIA4686+FeIIIA4658 and Fell+Felll A5262,5270 blends and line groups. Triple gaussian 
was used to deblend Ha with [Nil] A6548, 6583 doublet (for the components of the latter fixed 
wavelength difference and flux ratio F(A6583/A6548) = 3 were adopted). Exact (accuracy 
~ 0.1 A or better) wavelengths were taken from Coluzzi (1996). 

Table 2 lists the parameters of all the detected emission lines in the extracted spectra. 
All the lines redward of 5700 A are from the MPFS data, most of the bluer lines have fluxes 
measured with SCORPIO. Radial velocities from the panoramic data are given when possible. 
Number of lines such as [N I]A5200 and [S II]A4068 are close doublets or blends with unknown 
component ratios and therefore may have biased radial velocities by several tens of km s _1 . 
Errors in table 2 give the uncertainties of gaussian parameters. Signal-to- noise ratio is S/N ~ 
3 for the faintest lines identified like Fe IIIA4936. All the fluxes in the second column are 
normalized over the integral H/? flux in the corresponding spectrum, and unreddened fluxes 
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over the unreddened H/3 flux value. Total flux from the MPFS data is F(H/3) = (4.66 ±0.13) x 
ig-15 er g cm -2 g -i Unreddened values are calculated using the reddening curve by Cardelli et 
al. (1998) with Ry = 3.1 and the interstellar extinction value Ay = l m .54 calculated in section 
3.2. Unreddened H/3 flux is F (unreddened) (H/3) = (2.50 ±0.07) x 10" 14 erg cm" 2 s" 1 . 

Single strong lines show non-gaussian profiles possibly broadened by kinematical effects 
of the order 200 km s -1 . We however lack spectral resolution to deblend the broad and narrow 
components detected by Dunne et al. (2000). Some of the line identifications are debatable, 
mostly for weaker lines. However, we definitely see a number of high-excitation emissions like 
Pickering He IIA5412 and two putative [Ar IV] lines. Lines of moderately high excitation like 
He I, Fe III and S III dominate in number among the weaker lines. Some of the emissions are 
forbidden and therefore of nebular origin. FWHM of the He IIA4686 emission does not differ 
significantly from these of the brightest nebular lines like [O III]A4959,5007 hence He IIA4686 
is likely to be a nebular emission as well. Stellar or accretion disc emission is expected to be 
broadened by about 1000 km s" 1 (Hamann et al. 1995). Wolf-Rayet N III-N IV and CIII/CIV 
blends at about 4640 A are definitely absent (at least fainter than the Fe IIIA4658 emission) 
also supporting nebular origin of the He IIA4686 emission. 

3.2. Characteristic Line Ratios 

Ha/H/3 flux ratio is commonly used as a probe for interstellar absorption because of 
its very little dependence on all the plasma parameters. The ratio is nearly constant and 
equal to 3 with ~10% accuracy in the 5000 — 20 000 K temperature and 1 — 10 4 cm -3 number 
density range for photoionized plasma. We used Cardelli et al. (1998) extinction curve with 
R = 3.1. The integral MPFS spectrum gives Ay = l m .34 that is slightly greater than the Galactic 
absorption Ay al = l m . 14 in the X-ray source direction (Schlegel et al. 1998). Ha/H/3 ratio might 
be higher than 3 (Dopita & Sutherland 2003) in case of collisional excitation. The effect may be 
important, for example, for shock waves of moderate velocities Vs ^5 100 km s~ 1 (Shull & McKee 
1979), but its impact is questionable and virtually unstudied. Using the electron temperature 
and density estimates calculated below this section (T ~ 15 000 K, n e ~ 600 cm -3 ) one arrives 
to a slightly lower Ha/H/3 ratio ~ 2.8 that corresponds to Ay ~ 1™54 rather than the value 
stated above. 

Previous spectral studies of MF16 by Blair & Fesen (1994) and Blair et al. (2001) 
report higher values of interstellar absorption: E(B — V) = 0.52 (implying Ay ~ 1™6) and 
E(B — V)= 0.65 (Ay ~ 1™8), correspondingly. Possible explanation is in patchy extinction 
making the Ha/H(3 flux ratio variable within the extent of the nebula. This is consistent 
with the relatively low intrinsic absorption (intrinsic colour excess E(B — V) < 0™ 1, implying 
Ay < m .3) estimated for nearby stars a-c (Blair et al. 2001) and high absorbing column obtained 
for the X-ray source by Roberts & Colbert (2003) Nh ~4 x 10 21 cm -2 corresponding to Ay ~ l m .8 
(Gorenstein 1975). 
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Fig. 2. MF16 integral spectra extracted from the SCORPIO long-slit data. 
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In order to probe the physical conditions in the nebula we used several char- 
acteristic emission line ratios. Characteristic line calibrations were made according 
to Osterbrock & Ferland (2006) and STSDAS on-line NEBULAR.temden resource at 
http:/ / stsdas. stsci. edu/ nebular/ temden. html. 

[S II] 6717 / 6731 intensity ratio is usually used for electron density estimates. Integral 
line flux ratio (1.05 ± 0.03) indicates a rather high value, n e = 570 ± 60cm -3 if one assumes 
T= 15000K (see below). The electron temperature is usually estimated using two characteristic 
line ratios: [O III] (A4959 + A5007) /A4363 and [N II] (A6548 + A6583) /A5755. For different 
ions we estimate electron temperature as: T([0 III]) = 17700 ± 1200K, T([N II]) = 15600 ± 
2000K. Somewhat hotter temperatures inferred from [O III] lines are quite reasonable, because 
[O III] has higher ionization potential and emits in regions with higher temperature. An 
estimate for the coldest regions of the nebula may be made using [S II] (A4068 + A4076) / (A6717+ 
A6731) flux ratio. Using temden we arrive to a temperature estimate T([S II]) = 8600±500K. 
Large wavelength difference makes this value less defined due to uncertainty in the interstellar 
absorption value, T([S II]) = 9 000 ±1 000 K. Here, absorption uncertainty ~ 0™2 results in 
about 12% flux ratio uncertainty for the temperature-sensitive doublets of [S II] making total 
flux ratio uncertainty about two times higher. Accurate temperature determination is limited 
by the variable absorption effects mentioned above. 

Total luminosity of the nebula in the H/3 emission line may be used to estimate the 
emission measure and the total mass of the emitting gas. If one assumes constant density and 
temperature of the order (1 — 2) x 10 4 K for all the ionized material the mass of the nebula may 
be expressed as follows: 

yr ^ mEM ^ mL{H/3) /,n 

iV1 — n e E(Hf3) — aeff(H/3)n e E(H/3)> ^ 

where m is the mean particle mass and a e ^ is the effective (Case B) recombination 
coefficient. L(Hj3) = (9.0 ± 0.2) x 10 37 erg s _1 is the object luminosity is H/3. Using the 
estimates for n e = 600cm" 2 and T e = 15000K one yields for the case of MF16 M = 870 ± 50 M Q , 
about two times higher than the value obtained by Blair et al. (2001) because of higher H/3 
luminosity inferred. That is similar to the mass of the shovelled material for a spherical nebula 
with R — 13 pc and original interstellar hydrogen density no — 6cm -3 . 

4. Ionization sources 

4-1. Shock Excitation 

Most of the observational properties of MF16 such as two-component emission lines 
(Dunne et al. 2000), morphology (Blair et al. 2001) and high [S II]AA6717,6731 / Ha flux ratio 
favour shock excitation. Dopita & Sutherland (1996) provide a useful formula for the total flux 
emerging from a unit radiative shock front surface area in H/3 (both for the post-shock flow 
and for precursor): 
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1 Hf3 Fh/3, shock ± Fj] f3 , precursor 

a f 1/ \ 2.41 

7.44 x 10~ 6 ' v 
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100 km s / , . 

+9-86 x 10- 6 (too^hft) )x 
x 1 n ° o erg cm -2 s" 1 

lcm _,) ° 

Here no is the pre-shock hydrogen density and V is the shock velocity in 100 km s _1 
units. The estimates are supposed to be valid for shock velocities of the order 100 — 1000 km s _1 
and densities usual for the ISM and nebulae (0.1 — 10 4 cm -3 ). Precursor emission (the second 
term in brackets) is present only for shock velocities Vs ^ 150 km s _1 because slower shocks 
are unable to produce an ionization front moving fast enough (Dopita Sz Sutherland 2003). 

Integrating (2) over shock fronts gives total H/3 luminosities that may be compared with 
the observed luminosity of the emission-line object. The integral unreddened H/3 luminosity 
following from the 3D data analysis (see section 3.1) is L(Hj3) = (9.0 ± 0.2) x 10 37 erg s" 1 or 
(2.50 ±0.07) x 10~ 14 erg cm -2 s _1 in flux units. Formula (2) may be rewritten as: 

F{H(3)= (2.0 ( 1 5oT^1Ft) 2 ' 41 + 

+ 2 4f Y. ) 228 \ no Q2 y (3) 

+ V250 km s- 1 ) J 5cm-3 ° X V ' 

x 10~ 15 erg cm -2 s" 1 

9 is the angular diameter of the nebula in arcseconds. The estimate (3) predicts B./3 flux several 
times less than the observed value. Shock waves are short in producing enough flux in B./3 unless 
the initial density no > 10 cm -3 (that is unlikely because the mass of the shocked material should 
be much higher than the mass of the emitting gas) or the shock velocity > 300 km s -1 (that 
is excluded by kinematical data). Besides this, observed intensities of [O III] and [N II] lines 
require shock with Vs ~ 300 — 400 km s~ x (Dopita & Sutherland 1996; Evans et al. 1999) that 
is significantly higher than what one may allow having the kinematical limits mentioned above. 
He IIA4686 emission line is much brighter with respect to H/3 than a shock wave with precursor 
with Vs < 1000 km s _1 can provide (Evans et al. 1999). We conclude that observed expansion 
velocities are too low to explain the emission line luminosities observed as well as the line ratios. 

4-2. Photoionization 

An evident solution is to consider a bright photoionizing source that can increase the 
total luminosity in H/3 without changing the line widths. One may suggest the X-ray radiation 
of the central source responsible for powering the nebula. In order to check the ability of the 
X-ray source to power the nebula we set the radius of the nebula equal to 13pc (according to 
Blair et al. (2001)) and used the best-fit MCD + power-law model from Roberts & Colbert 
(2003). MCD (MultiColour Disc blackbody) was calculated according to Mitsuda et al. (1984). 
Soft MCD component was extrapolated to the EUV region, but we truncated the power-law at 
0.5 keV where its intensity becomes comparable with that of the MCD component. 
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The observed X-ray radiation provides an ionization parameter logarithm \gU ~ —5.2 
2 . However, observed [O III]A5007/H/3 ~ 7 and HeIIA4686/H/3 ~ 0.2 ratios point to a rather 
high ionization parameter logarithm value IgU > —3 (Evans et al. 1999). Besides this, as we 
see below, the X-ray source itself is unable to produce enough photoionizing quanta. This can 
be done by introducing a powerful extreme ultraviolet (EUV) source (Ljjy ~ 10 erg s ). 

He IIA4686 emission line is much brighter with respect to H/3 than a shock wave with 
reasonable parameters may explain. Total luminosity in the line is Lh c iia4686 = (2.0 ± 0.2) x 
10 37 erg s~ x according to the MPFS data. Because of this and the high excitation potential of 
the corresponding transition, we consider He IIA4686 a recombination line. One can estimate 
the number of ionizing quanta required and the luminosity of the photoionizing source following 
Osterbrock & Ferland (2006) (Zanstra method): 

Q(He + ) > £ (A 4 686) x 

> V//(^A4686) L(HeIIX468Q), 
as and a e ^ (Hel I A4686) here are the total recombination rate for the Case B and the effective 
recombination rate for the HeIIA4686 line (Osterbrock & Ferland 2006), correspondingly. Their 
ratio depends weakly on plasma parameters, varying by a factor of 2 in the plasma temperature 
range (1 — 30) x 10 3 K. For particle density of the order 100 — 1000 cm -3 and temperature T ~ 
(1.5 — 2) x 10 4 K the recombination coefficient ratio as /<y e ^ (Hell X4686) ~4.7. Corresponding 
luminosity of the He+-ionizing EUV source is L(A < 228 A) > 95L(HeIIA4686) ~2x 10 39 erg s" 1 . 
Note that all the ionizing luminosity must be concentrated in the wavelength range 20 — 200 A 
(at A ^ 20 A heavier elements dominate the extinction). The luminosity is calculated for 
A = 228 A and actually is a lower estimate because of higher photon energies and quanta leakage. 
Though bright HeIIA4686 emission seems to be quite usual for ULX counterparts (Abolmasov 
et al. 2007), for MF16 the Hell luminosity is extremely large. The EUV luminosity inferred is 
comparable with the X-ray luminosity. 

Similar estimates may be made for other recombination lines such as H/3 and HeIA4471. 
The Case B numbers of He + , He and iJ°-ionizing quanta derived from HeIIA4686, HeIA4471 
and H/3, correspondingly, are as follows: 

Q(He-) = (2.2 ± 0.2) x 10* ^^4686) (5) 
VV ' V ; 2.0 x 10 37 erg s -1 V ; 

/ ru / 4Q L(HeIXU71) i 

Q (He ) = 1.8 ± 0.2 x 10 49 — s" 1 6 

4.0 x 10 36 erg s" 1 

Q(H°) = (1.90 ± 0.05) x 10 50 T s- 1 (7) 

vv ; v ; 9.0 x 10 37 erg s" 1 V ; 
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Here and below we define ionization parameter as (following Evans et al. 1999): 

i r+°° p 

U= / £-d{hv), 

cnn J13.60V hv 

where nn is the total (ionized + neutral + molecular) hydrogen density in the gas. 
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The ionizing quanta numbers listed above correspond to ionizing source luminosities greater 
than 10 39 erg s" 1 . Luminosity required to explain the Hj3 luminosity of MF16 is L(X < 912 A) > 
4 x 10 39 erg s^ 1 . 

4-3. Photoionization grid 

We have computed a grid of Cloudy version 07.02.00 (Ferland et al. 1998) photoionization 
models in order to fit the spectrum of MF16 neglecting shock waves. We assumed all the plasma 
situated at 13 pc from the point source in the center, forming a hollow envelope with hydrogen 
density rtu = 500 cm -3 . We considered photoionization with an isotropic source with a two- 
component spectrum: fixed X-ray component observed by Chandra (and extrapolated in the 
EUV region in the same way it was done in the previous section) and a blackbody source with 
the temperature spanning the range from 10 3 to 10 6 K and integral flux densities from 0.01 to 
100 erg cm -2 s -1 (corresponding to the luminosity range 10 37 — 10 41 erg s -1 ). 

8 lines with different ionization potentials were used for fitting. Measured emission line 
fluxes were fitted assuming the distance D = 5.5 pc. Two interstellar extinction values were 
applied, Ay = l m . 34 and 1™ 54, the latter resulting in much better fit supporting the idea that 
the flux ratio Ha / H/3 ~ 2.8 rather than 3. Fitting results and partial x 2 values are given in 
table 3. 

The best-fit parameters for MF16 are: lgT(i^) = 5. 15 ±0.10, F = 0.6 ±0.1 erg cm -2 s" 1 
resulting in a normalized x 2 — 25 (for 6 degrees of freedom). The luminosity of the corresponding 
UV source is Ljjy ~ 1.2 x 10 40 erg s _1 (mainly concentrated at 200 — 400 A). The high values of 
the reduced x 2 are a direct result of the high S/N in our data. Note that all the line fluxes are 
predicted with the accuracy better than ~ 20% quite reasonable for our analysis. All the line 
ratios are affected by plasma inhomogeneities, abundance and depletion effects. In principle 
a more complicated model may be applied to our data. Hotter lines like He IIA5412 and 
[O III]A4363 are enhanced in the observed spectrum with respect to the best-fit photoionization 
model predictions indicating an additional source of heating. Model [O III] temperature is 
about 10 4 K instead of 1.7 x 10 4 K measured from the [O III] characteristic line ratio. 

Our models also overestimate the intensities of Fe III lines. This may be attributed to 
partial depletion of the element into dust taking place in the interstellar medium. Iron depletion 
values are usually of the order 10 -2 for undisturbed dense ISM and about 0.1 for HII- regions 
(Rodriguez 2002) indicating that dust is effectively destroyed in MF16. Shock waves are likely 
to destroy significant part of dust grains in the ISM but some percentage of grains survive for 
Vs ^ 500 km s~ x (Dopita & Sutherland 2003; Jones et al. 1994). The rich spectrum of Fe III 
indicates that the emitting gas was processed by shock waves showing however depletion by a 
factor of 0.3 - 0.5. 

We conclude that the observed emission-line spectrum may be reproduced without any 
significant abundance changes and without shock waves. However, the effect of the latter may 
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still be important (in heating the medium and destroying dust). 
5. Discussion 

5.1. The Nature of the EUV Source 

Cool multicolour disc (MCD) + power law models are often used to fit ULX spectra 
and are considered an argument for the IMBH model (Colbert & Miller 2005). Hard power 
law component usually dominates the X-ray spectrum but it is expected to be truncated at 
energies close to the inner disc temperature. For NGC6946 ULX-1 about 20% of the standard 
Chandra X-rays may be attributed to the thermal MCD component with the inner temperature 
T in = 0.15 keV (Roberts & Colbert 2003). Bolometric luminosity of a standard disc with such 
a temperature should be about 5 times higher than its luminosity in the standard Chandra 
X-ray range. Though most of the disc luminosity is emitted in the EUV range the number of 
photons is insufficient to explain the Balmer lines. It may be checked that a conventional MCD 
with the parameters measured by Chandra provides only Qh — 10 49 s _1 - at least an order of 
magnitude less than is needed to explain the H/3 luminosity. Luminosity of the He IIA4686 line 
is however only about 2 times higher than the MCD component can provide. 

The central optical point source (star d in figure 1) is much brighter than a possible low- 
energy tail of a standard accretion disc. According to Blair et al. (2001), star d has V = 22 m . 64 
and B — V = n !46. For a distance of 5.5Mpc and Ay = 1 II !54 the central optical point-like 
source will have My = —7™ 6 (Ly ~ 3 x 10 38 erg s -1 ) instead of My = —2™ 7 given by the 
best-fit MCD for Chandra data with an infinite outer radius. The optical excess, however, 
may be attributed to the donor star. The optical object itself is relatively cool: (B — U) ~ 
— n ?05 if one assumes colour excess following from the Ay estimates due to Ha/H/3 ratio 
E(B — V) ~ Ay /3.1 ~ 0.5. That corresponds to a temperature ~ 10 A K. Exponentially small 
part (exp{— Ry/T} ~4x 10 -5 ) of the photons emitted by the star is able to ionize hydrogen 
and helium, therefore the contribution of the donor star to the ionization balance in MF16 is 
probably negligible. About l m additional intrinsic absorption is needed to make the optical 
source a hot blackbody with B — V ~ — 0.3 m . 

In figure 3 we reconstruct the spectral energy distribution (SED) of the ULX throughout 
the UV spectral range (from X-rays to the optical). Our estimates of ionizing luminosities 
put together with the X-ray best-fit model spectrum (Roberts & Colbert 2003) result in a 
source with an approximately flat [yL v ~ const) spectrum from 0.01 to 10 keV. Ionizing 
flux estimates (5-7) are indeed lower limits due to higher energies of the absorbed quanta and 
possible quanta leakage. The estimated EUV luminosity can not be explained by an MCD with 
Ti n ~ 0.1 —0.2 keV. In figure 3 we show MCD spectra for different black hole masses (10, 100, 
1000 and 10000 M from bottom left to top right) accreting at 1% of the critical accretion rate. 

For a standard disc (Shakura & Sunyaev 1973) luminosity and temperature scale with 



13 



the black hole mass M and dimensionless accretion rate m as Lot Mm, Tj n oc M _1//4 m 1//4 , 
correspondingly. Mass and accretion rate estimates are therefore much more sensitive to the 
SED shape than to luminosity: M oc L 1 ' 2 T j ^ 2 , m oc L x ' 2 Tf n . The black hole mass in the 
framework of the IMBH hypothesis should be ~ (1 — 3) x 1O 4 M . 

Our best-fit model parameters (section 4.3) suggest an object with radius ~ 10 R Q and 
temperature ~ 10 5 K too high even for a conventional Wolf-Rayet star. Hot WC and WO 
stars may have comparable effective temperatures, but their luminosities and ionizing quanta 
production rates are at least an order of magnitude lower (Crowther et al. 2006). A possible 
interpretation is emission from the optically thick wind of a supercritically accreting stellar 
mass black hole that should be a bright UV/EUV source of appropriate luminosity (Poutanen 
et al. 2007). The X-ray component may be attributed to the funnel radiation in that case. 
Poutanen et al. (2007) estimate the outer photosphere temperature of a supercritical accretion 
disc wind as: 

T ph ~ 0.8m- 1/4 m" 3/4 keV (8) 

where m is the black hole mass in solar units, rh is the dimensionless mass accretion rate at 
the infinity. From the equation (8) the dimensionless mass accretion rate may be estimated as 
m ~ 200m 1 T 5 4 ^ 3 , where mi is the black hole mass in 10 M Q units, and T5 is the temperature 
of the photosphere of the wind in 10 5 K. The accretion (or, more strictly, mass ejection) rate is 
found to be about an order of magnitude less than expected for SS433 (Fabrika 2004) yet still 
highly supercritical. 

5.2. ULXs m the UV 

Let us assume the SED of the central object smooth throughout the UV range and 
interpolate betweeen the optical and EUV flux estimates with a power-law (see figure 3). The 
predicted intrinsic (unabsorbed) flux at 1000 A is F\ ~ 10~ 15 erg cm -2 s _1 A -1 . At shorter 
wavelengths radiation is highly absorbed by neutral gas. Absorption by dust at A ~ 1000 — 
2000 A may be estimated as about (3 — 4) x Ay ~ 5 m (Cardelli et al. 1998), so one should expect 
at 1000 A Fa ~ 10~ 17 erg cm~ 2 s~ x A~ x , that corresponds to uiab ~ 25 m (we use tuab definition 
from Oke & Gunn 1983). 

Predicted UV flux is quite reachable for GALEX pointing observations but too faint 
to be detected in the All-sky and Medium-sky imaging surveys (Martin et al. 2005) mainly 
because of high interstellar absorption in UV spectral range. Less absorbed ULXs might be 
quite achievable targets if they are as bright as NGC6946 ULX-1 in the UV. In figure 3 we 
compare the vF v values with flux limits (10 5 s, S/N = 5) for GALEX photometry in FUV (1400- 
1700 A) and NUV (2000-2700 A) bands. Signal-to-noise estimates were made with the GALEX 
Exposure Time Calculator (http://kaweah.caltech.edu:8000/ExpCalc.tcl). 

Capabilities of the HST instrumentation may be used as well to study ULX sources, at 
least those less absorbed than NGC6946 ULX-1. Liu et al. (2005) studied the UV counterpart 
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of NGC5204 X-l (Ul object) with STIS MAMA in the FUV range and classified the UV object 
as a BOIa star with some oddities like strong NVA1240 emission that may be attributed to an 
ambient HII region, irradiated accretion disc or corona. Some resonance lines like SHVA1400 
and CIVA1550 that show P Cygni-type profiles in OB supergiant spectra are missing (or present 
only as weak absorptions) in the spectra. Liu et al. (2005) explain this effect by Roche lobe 
overflow by the B supergiant donor star in a binary system. This argumentation seems to be 
questionable because the winds of early B-stars have terminal velocities about 1000 km s _1 , an 
order of magnitude higher then virial velocities characteristic for HMXBs. There is however 
an alternative explanation for the unusual absorption spectrum of Ul: P Cygni profiles of UV 
lines formed in hot rarefied atmospheres of OB Ia/b supergiants transform into absorptions 
with increasing mass loss (Hutchings 1976). Mass loss rate needed for these effects to become 
important is of the order 10 -4 M yr _1 . It is also possible that the OB-supergiant spectrum 
originates not from the donor star but from the accretion disc wind atmosphere, that is probable 
if the mass accretion rate is close to that observed in SS433. 

Monochromatic flux at 2200 Areported by Liu et al. (2005) for Ul is F\ ~ 3.7 x 
10~ 17 erg cm -2 s -1 A. The best-fit absorbing column measured for the X-ray source implies 
E(B -V) = m 28. Galactic absorption is negligibly small (E(B - V) = m 013 according to 
Schlegel et al. (1998)). Intrinsic UV luminosity of NGC5204 X-l calculated using these val- 
ues appears to be vL v ~ 10 38 erg s -1 not taking into account absorption or ~ 10 39 erg s" 1 if 
reddening calculated from X-ray fitting is used. 

Unfortunately, STIS stopped science operation in 2004 and its potential for studying 
ULXs in the UV was lost. ACS SBC and HRC cameras require exposures > 10 5 s to obtain 
S/N> 1 (exposure time calculator at http://apt.stsci.edu/webetc/acs/acsJmg-etc.jsp was used) 
in the case of NGC6946 ULX-1 but may be used to study less absorbed bright UV counterparts. 
Among the most perspective ULXs are NGC5204 X-l, NGC4559 X-7 (Soria et al. 2005), HoII X- 
1 (Lehmann et al. 2005) and HoIX X-l (Miller 1995) but we do not have any direct evidence 
that these objects are as bright as NGC6946 ULX-1 (line luminosities are generally lower). 

The two most popular models of ULXs (IMBHs and supercritical accretion discs) both 
predict strong UV/EUV emission. IMBHs are expected to be bright UV sources due to the disc 
temperature dependence on the mass of the accretor (T in cx M _1//4 ). The greater is the mass of 
the black hole the higher is the UV flux, approximately as Foe at a given wavelength for 
fixed dimensionless accretion rate (Shakura & Sunyaev 1973; Mitsuda et al. 1984). Standard 
accretion discs around IMBHs are expected to be truncated only at very low temperatures, 
close to 1 eV (if their outer radii are due to tidal truncation; see, for example, discussion in 
Abolmasov et al. (2007)), and must have therefore a distinguished spectral slope characteristic 
for an MCD with an infinite outer radius, F v oc z/ 1//3 (Shakura & Sunyaev 1973). 

Advective super-Eddington discs are often fitted by slim disc models (Abramowicz et 
al. 1988) predicting roughly flat SED in the X-ray band. On the other hand, supercritical 
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accretion discs are expected to have strong optically thick winds (Shakura & Sunyaev 1973) 
with photospheres having blackbody-like spectrum peaking somewhere in the EUV/UV region 
(Poutanen et al. 2007). The exact structure of the accretion flow is irrelevant for the EUV/UV 
SED because most of the EUV radiation comes from the wind photosphere. In the GALEX 
range the spectrum is expected to have a Ray leigh- Jeans shape F v oc v 2 . 

We conclude that detecting ULXs in the GALEX UV range and measuring the spectral 
slope may help to distinguish between these two models of ULXs (IMBH or supercritical accre- 
tion). At least for some subsample of ULXs (those having indications for a hard photoionizing 
continuum) different models predict different spectral slopes in the UV range (see figure 3). 

The nebula may contribute to the UV emission from the source. Cloudy model predicts 
total flux in the GALEX range about an order of magnitude less than that expected from the 
ULX itself. Best-fit model predicts L(Lya) ~ 2QL(H(3) ~ 2 x 10 39 erg s _1 but most of the 
La-quanta are likely to be absorbed by neutral interstellar medium. Expected contribution 
from the nebula is shown in figure 3. Higher fluxes may appear if some additional amount of 
gas with at least comparable emission measure and much higher temperature is present (much 
cooler gas does not contribute to the continuum in the UV range). But the optical spectra 
contain no signatures of high-temperature plasma such as coronal lines of very high-ionization 
element species. Therefore nebular emission is unlikely to give a significant contribution to the 
GALEX UV in the particular case of MF16. 

5.3. MF16 as a Jet-Blown Shell 

From the high spatial resolution HST images (Blair et al. 2001) it can be seen that the 
nebula is elongated with an axis ratio no less than ~ 1.5. Either an underlying density gradient 
or highly anisotropic energy injection in the ISM is needed to explain the observed morphology. 

If a shell evolves in a strong density gradient, one side expands with a much higher 
velocity making the nebula much more extended in one direction then in the other. Central 
position of the optical star d is difficult to explain in this scope. Density gradient probably plays 
significant role only in making one side of MF16 much brighter. This is expected (Maciejewski 
& Cox 1999) even for very shallow density gradients. 

In all the self-similar solutions usually applied to expanding shells (Lozinaskaya 1992) 
and wind-blown bubbles (Castor et al. 1975), energy and mechanical luminosity appear with 
very low exponents. In particular, for the Sedov solution for an adiabatic-stage SNR the radius 
of the shell depends as R oc E 1 ^ on the energy input. Therefore to achieve 50% variation in 
the radius of the nebula, energy injected per unit solid angle (or the initial ISM density) must 
change no less than by an order of magnitude. However, if the central machine is a source of 
well-collimated wind or jets, the shape is naturally explained. 

Jet activity is unlikely to be present in standard discs (expected in the framework of 
IMBH model) and may be considered an argument for supercritical accretion. We do see a pair 
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Fig. 3. Reconstructed SED of NGC6946 ULX-1 from X-rays to the optical. Part of the best-fit X-ray 
spectrum (Roberts & Colbert 2003) is shown by a thick solid line at shorter wavelegths. Dashed 
curve shows the behavior of the Chandra MCD component (outer radius is set to infinity) at longer 
wavelengths. Ionizing flux estimates for Hell, Hel and HI are shown by diamonds with error bars. 
Asterisk corresponds to the optical source d (Blair et al. 2001) in the V band corrected for inter- 
stellar absorption Ay = 1™54. MCD SEDs are shown (thin solid lines) for m = 0.01 and masses 10, 
100, 1000 and 10 4 M a . The best-fit black body obtained during our Cloudy modelling (see sec- 
tion 4.3) is shown by a dotted line. GALEX photometric limits (FUV and NUV bands) are given 
for S/N=5, exposure time 10 5 s. Gray line shows the best-fit Cloudy model spectrum of the nebula. 

of powerful jets and an elongated nebula W50 in the case of SS433. However, if NGC6946 ULX- 
1 is a supercritical accretor viewed face-on (at inclination i < 20°) the direction of jets must be 
close to the line of sight. In that case the actual length of MF16 should be several times (about 
sin _1 z) greater, approaching the size of W50. 

NGC6946 is a seen nearly face-on: according to Gordon (1968), its inclination is in the 
range 16 — 27°. If the actual size of MF16 in the direction perpendicular to the disc plane of 
the host galaxy is of the order 100 pc as in the case of W50 then the nebula is likely to evolve 
in an ISM with a considerable density gradient. The thickness of the gaseous disc in NGC6946 
probably does not differ too much from the Galactic neutral hydrogen disc thickness ~ 300 pc 
(Lockman & Gehman 1991). The expected effects are similar to those observed in the case of 
W50 that evolves in a strong ISM density gradient (Lockman et al. 2007). It is quite probable 
that the asymmetry of MF16 seen in HST imagesoriginates from the interaction of one of the 
jets with denser material closer to the galactic plane. 
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6. Conclusions 



In our optical spectra of MF16 we detect more than 30 emission lines including several 
high-excitation lines of He II and [Ar IV] and a rich spectrum of moderately high excitation 
Fe III lines. He II lines are narrow (broadened by < 300 km s" 1 ). We also do not detect any 
Wolf-Rayet features. Therefore we suggest that the He IIA4686 line is of nebular origin. 

Large intensities of iron lines indicate moderately low depletion of the element into dust 
(30 — 50% in the gas phase), probably due to dust destruction in shock waves. Destruction is 
probably incomplete because Cloudy models overestimate the intensities of Fe III lines. 

We find the electron density in the nebula n e = 570 ±60 cm -3 ([S II] A671 7,6731), electron 
temperatures for different ions are T([0 III]) = 17700±1200K, T([N II]) = 15600±2000K and 
T([S II]) = 9000±1000K indicating the presence of regions with different electron temperatures. 
Total hydrogen emitting gas mass is M ~ 900 M . Interstellar absorption is Ay ~ 1"?34 
in traditional Ha/H/3=3 assumption. More realistic value Ha/TT3=2.8 results in a higher 
extinction value Ay = 1™ 54. 

The observed line luminosities and diagnostic line ratios appear to be inconsistent with 
excitation and ionization by shock waves, therefore we suggest an EUV source responsible 
for powering the nebula. Photoionization modelling with Cloudy as well as Zanstra estimates 
suggest the central source must be ultraluminous not only in X-rays but also in the UV/EUV 
range emitting about 10 40 erg s _1 in the spectral range 100 — 1000 A. Using the observed 
X-ray spectrum, ionizing flux estimates from He II, He I and H lines and the optical point-like 
counterpart (star d) we reconstruct the SED of the central object from X-rays to the optical. 
The derived spectrum is roughly flat, vL v ~ const. 

Both most popular models of ULXs (IMBHs and supercritical accretion discs) predict 
high UV/EUV luminosities. However these two models predict different spectral slopes in the 
UV region. We conclude that measuring the spectral slope will help to distinguish between the 
two models as well as to determine the parameters of the successfull model such as the black 
hole mass in the case of IMBH or accretion rate in the case of a supercritical accretion. 

Actual excitation and ionization conditions in MF16 may be much more complicated. 
Though introducing a bright EUV source is a possible way to explain the spectrum, we suggest 
that other explanations such as high pre-shock density and multiple shocks with different ve- 
locities are not completely excluded. Observations with higher angular and spectral resolution 
(but also higher S/N ratio than that used by Dunne et al. (2000)) may provide additional 
information about the enigmatic nebula MF16. 
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Table 1. Observational logs 



Spectrograph MPFS 

Date 2005/17/01 

Total exposure (s) 5829 

Spectral range ( A) 4000-7000 

Spectral resolution (A) 6 

Seeing, arcsec 1.7 



SCORPIO (long-slit mode) 
2005/10/06 
5400 
3900-5700 
5 
l.G 
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Table 2. MF16 emission line parameters. Line intensities in H/3 units (the uncertainties do not account for the the H/3 flux 
uncertainty), radial velocities and FWHMs. 



line V^kms' 1 FWHM,A 

(unreddened) 



[S7/]A4068 


0.200±0.013 


0.276±0.018 


-100±20 


8.3±0.7 


HS 


0.228±0.013 


0.311±0.018 


-66±11 


6.6±0.4 


[Tim] A4161? 


0.018±0.008 


0.024±0.010 


0±40 


3.8±1.2 


S777A4284? 


0.074±0.012 


0.094±0.015 


-50±40 


11.2±1.3 


#7 


0.408±0.013 


0.510±0.016 


-47±14 


6.2 ±0.2 


[OJJJ]A4363 


0.146±0.016 


0.181±0.019 


-73±14 


7.1±0.5 


[Fe77]A4414? 


0.039±0.010 


0.047±0.012 


-10±50 


8.6±1.6 


HeIXU71 


0.038±0.006 


0.045±0.007 


-56±20 


6.0±0.7 


[Fe/U]A4658 


0.066±0.011 


0.072±0.012 


10±30 


9.3±1.2 


He7/A4686 


0.205±0.009 


0.222±0.010 


-5±18 


6.2±0.2 


[FeIII]XA701 


0.017±0.008 


0.018±0.008 


-35±50 


5.1±1.7 


HeJA4713 (+[Ar/y]A4711?) 


0.021±0.008 


0.021±0.008 


-70±40 


5.3±1.6 


[ylr/V]A4740 


0.009±0.004 


0.009±0.004 


-50±30 


3±1 


H/3 


1.000±0.017 


1.000±0.017 


-28±4 


6.2±0.2 


[FeIII\X4S81 


0.018±0.008 


0.018±0.008 


-86±20 


4.8±1.0 


ifeIA4922 


0.011±0.005 


0.011±0.005 


-10±30 


3.9±1.2 


[Fe///]A4936 


0.008±0.005 


0.008±0.005 


-60±50 


4. Oil. 5 


[OIII] A4959 


2.337±0.050 


2.24±0.05 


-23±3 


6.2±0.1 


[FeIII]X4986 


0.044±0.010 


0.042±0.009 


50±100 


4.0±1.5 


[OJJJ]A5007 


6.94±0.14 


6.50±0.13 


-17.4±2.5 


6.2±0.1 


[Fe/II]A5033? 


0.013±0.005 


0.012±0.005 


50±30 


4.1±1.4 


[Fe//]A5158? 


0.073±0.005 


0.064±0.004 


11±9 


7.8±0.4 


[iV/]A5199 


0.152±0.005 


0.131±0.004 


-34±5 


7.1±2.0 


[Fe77]A5262? 


0.072±0.014 


0.061±0.012 


140±90 


18±2 


[FeJJ] + Fell 1X5270 


0.018±0.008 


0.015±0.007 


60±30 


5.1±1.7 


Fe/77A5300? 


0.011±0.005 


0.009±0.004 


100±50 


5±2 


ffeIIA5411 


0.021±0.005 


0.017±0.004 


-10±20 


5.6±1.0 


[7V//]A5755 


0.118±0.015 


0.084±0.011 


-70±30 


16.0±1.5 


He/A5876 


0.114±0.024 


0.078±0.016 


-9±20 


7.0±1.1 


[OJ]A6300 


1.421±0.083 


0.89±0.05 


25±9 


10.1±0.4 


[OI] A6364 


0.480±0.029 


0.296±0.018 


12±9 


9.7±0.4 


[iVII] A6548 


1.385±0.027 


0.822±0.016 


-18.7±1.3 


7.6±0.1 


Ha 


4.728±0.087 


2.80±0.05 


-15.7±2.1 


7.8±0.1 


[iV//]A6583 


4.156±0.082 


2.45±0.05 


-18.6±2.2 


7.6±0.1 


He/A6678 


0.060±0.032 


0.035±0.018 


-18±9 


14.0±5.2 


[SII] A6717 


2.457±0.024 


1.412±0.014 


-37.5±1.1 


7.7±0.1 


[S7/]A6731 


2.348±0.025 


1.346±0.014 


-32.8±1.2 


8.2±0.1 
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Table 3. Fitting emission line fluxes with Cloudy models. Emission line fluxes are given in 10 erg 



line 


Fobs 


Fmodel 


Ax 2 


-^obs 


Fmodel 


A X 2 


Ay, 




1.34 






1.54 




Fc IIIA4658 


11±2 


30 


90 


18±3 


35 


32 


He IIA4686 


44±2 


36 


16 


56±2 


44 


3G 


R/3 


200±3 


205 


2 


249±4 


246 


0.56 


[0 III]A5007 


1310±30 


1089 


54 


1620±30 


1520 


11 


He IA5876 


16±3 


27 


13 


20±4 


32 


9 


[N IIJA6584 


526±10 


484 


15 


611±12 


527 


49 


[S II]A6717 


304±3 


334 


100 


352±3 


367 


25 


[S II]A6731 


290±3 


316 


75 


335±4 


348 


11 


F, erg cm -2 s _1 


u.oo_ 05 






o.6i° :l 




T, 10 5 K 


c nc+0.15 
°- uo -0.05 






c- -,+0.15 
°' i -0.05 




X 2 /DOF 




62 
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